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ABSTRACT 

NGC 55 ULXl is a bright Ultraluminous X-ray source located 1.78 Mpc away. We analysed 
a sample of 20 Swift observations, taken between 2013 April and August, and two Chandra 
observations taken in 2001 September and 2004 June. We found only marginal hints of a 
limited number of dips in the light curve, previously reported to occur in this source, although 
the uncertainties due to the low counting statistics of the data are large. The Chandra and Swift 
spectra showed clearly spectral variability which resembles those observed in other ULXs. 
We can account for this spectral variability in terms of changes in both the normalization 
and intrinsic column density of a two-components model consisting of a blackbody (for the 
soft component) and a multicolour accretion disc (for the hard component). We discuss the 
possibility that strong outflows ejected by the disc are in part responsible for such spectral 
changes. 

Key words: accretion, accretion discs - X-rays; binaries - X-rays: galaxies - X-rays: indi¬ 
viduals (NGC 55 ULXl, XMMU J001528.9-391319) 


1 INTRODUCTION 

The population of Ultraluminous X-ray sources (ULXs) has grown 
significantly after the launch of the XMM-Newton and Chandra 
satellites. These observatories made it possible to carry out the de¬ 
tailed spectroscopic study of intrinsically powerful but faint point¬ 
like X-ray sources in nearby galaxies (< 100 Mpc). They were 
dubbed Ultraluminous X-ray sources bec ause they have a very high 
X-ray luminosity, above 10^® erg s“^ l lFabbianolll989h . and ex¬ 
ceed (although not necessarily all the time) the Eddington limit 
for spherical accretion of pure hydrogen onto a 10 Mq black hole 
(BH). Their position is not coincident with the nucleus of their 
host galaxy, ruling out a low luminosity Active Galactic Nucleus 
(AGN). While a significant fraction of ULXs can be associated to 
background AGNs or young supemovae interacting with the cir- 
cumstellar matter, the majority of them are now thought to be a 
pecu liar class of accr eting sources, either BHs of stellar origin 
(e.g. iKing et ^ 1200 ih acc reting above Eddington or Intermedi- 
ate Mass BHs (IMBHs, e.g. lColbert & Mushotzkvl[l9^ accreting 
sub-Eddington. Becaus e of the special cond i tions/environment re¬ 
quired to form IMBHs jMadau & Reesl200ll : IPortegies Zwart et ^ 
1 20041 : Ivan der Marell l20 04ll. the latter sceimio is considered a vi¬ 

able one for HLX-1 dFarrell et ^ l2009l : [^illat et al.llMIh . that 
may however be interpreted also as the nucleus of a stripped dwarf 


galax y dSoria et alj 1201 ll : IWebb et ^ l20ld : iMapelli e t al.| |20lj, 
2017) , and possibly for a h andful of other sources (e.g. Sutton et^ 
2012 : |Pasham et al]|2014ll . On the other hand, for the great major¬ 


ity of ULXs an interpretation in terms of a stel lar mass BH or a BH 
of ste l lar origin seems mo re likely (see e.g. IZampieri & Robe^ 
1 20091 : iFeng & Sorial 1201 ll) . Recent dynamical measurements or 
constraints to the BH mass have been deriv e d and are in agree - 
ment with this conclusion dUiu et al.l 1201 3l : iMotch et alj l2014h . 
In addition, the observed X-ray spectral and variability proper¬ 
ties are consistent in most cases with an interpretation in terms 
of super-Eddington accretion and possibly beamed emission (e.g. 
iGladstone et al] l2009l : ISutton et alJ 120131) . as expected for light 
BHs, with the actual BH mass (5-80 M (7)) depending somewhat on 
the metallicity of the environment (e.g. Zampieri & Robert^l20091 : 


IMapelli et alj200lj : lBelczvnski et al.l2010h . Finally, it was very re¬ 

cently showm_tha^sometransient ULXs may even host Neutron 
Stars dBachetti et alj2014l) . 


One of the main diagnostics to study ULXs is the X-ray spec¬ 
trum. High quality ULX spectra are usually characterised by a 
soft component, well fit by a standard disc spectrum, and an ad¬ 
ditional high energy component with a turn-over at 3-5 keV, phe¬ 
nomenologically desc ribed in terms of Comptonization in an op- 
tically thick medium dStobbart et alj|200^ : IGladstone et al.ll2009l : 





















































2 F. Pintore, P. Esposito, L. Zampieri, S. Motta, A. Wolter 


IPintore & Zampierill2012h . Although this two-components model 
is only a phenomenological description, the peculiar properties 
of the ULX spectra has led to the conclusion that these sources 
are in a markedly different accretion regime from those observed 
in Galactic BH X-ray binary systems (XRBs) accreting at sub - 
Eddington rates ( [Gladstone et alJl200^ : lKaiava & Poutanenll200S}l . 
The combination of spectral shape, spectral evolution and variabil¬ 
ity suggest that the mas s accretion rate in the disc is at or above 
the Eddington limit l e.g. [Midd leton et al.ll201 iLISutton et al [20131 : 
(Bachetti et all[2013l: [Walton et al.[[2013[ : fpintore et alj [20141 . Un¬ 
der such conditions, strong, radiatively driven winds can originate 
from the inner regions of the disc a n d propagate outwards above 
its plane (e.g. ^utanen et alj [2007[ : [Ohsuga & Mineshi^ [20071 : 
[Ohsuga et alj200^ . These turbulent and optically thick winds may 
originate the soft spectral component observed at low energies 
and possibl y also the extrinsic v ariability detected at higher en¬ 
ergies (e.g. [Middleton et al.[|201 ih . while an advection-dominated 
disc, the emission of which may be comptonized, produces the 
high energy part of the spectrum (e.g. Middleton^tah[|2blT|. [20 1 j : 
[Sutton et al.[ [2OI3I : [Walton et al.[ 1^13 : Pintore et alJ 201 ^ . The 
structure and emission properties of the inner disc may be affected 
also by BH spin and radiative transfer effects taking place in th e 
disc atmosphere (e.g [Suleimanov et al.[ [2002[ : [Kawagu^ [2003h . 
This picture has been confirmed by observations ab o ve 10 keV 
recen t ly performed by NuSTAR ( [Bachetti et al.[ [20131 : [Rana et aP 
[2014 [Walton et al.[[2013[ . [20l4) . Clearly, the onset of strong out¬ 
flows would play an important role in modifying the structure of 
the accretion disc as they are expected to effectively remove energy 
and angular momentum from the system. 


Table 1. Summary of the Chandra and SWIFT observations of NGC 55 
ULXl. 


Chandra Obs. ID 

Start time 

Stop time 

Exp. 

Net counts 


(Terrestrial Time) 

(ks) 


2255 (Chip 10) 

2001-09-11 06:25:05 

2001-09-11 23:30:20 

59.4 

9509 

4744 (Chip 11) 

2004-06-29 01:48:01 

2004-06-29 05:17:04 

9.6 

3719 


Swift Obs.ID 


00032619001 

2013-04-10 03:24:17 

2013-04-10 09:50:56 

4.8 

301 

±17 

00032619002 

2013-04-17 11:22:00 

2013-04-17 16:30:56 

5.2 

299 

±18 

00032619003 

2013-04-24 13:35:09 

2013-04-24 23:19:54 

4.8 

178 

±14 

00032619004 

2013-05-01 00:48:38 

2013-05-01 05:54:54 

5.0 

191 

±14 

00032619005 

2013-05-08 05:45:04 

2013-05-08 23:54:57 

4.9 

258 

±16 

00032619006 

2013-05-15 18:48:28 

2013-05-15 22:27:55 

4.1 

167 

±13 

00032619007 

2013-05-22 15:47:03 

2013-05-22 22:49:55 

5.6 

182 

±14 

00032619008 

2013-05-29 19:32:23 

2013-05-29 19:39:54 

4.4 

160 

±10 

00032619009 

2013-06-02 11:33:13 

2013-06-02 15:04:55 

4.5 

360 

±19 

00032619010 

2013-06-05 00:18:17 

2013-06-05 05:26:57 

4.6 

483 

±22 

00032619011 

2013-06-12 18:28:52 

2013-06-12 23:25:56 

3.4 

168 

±13 

00032619012 

2013-06-19 18:28:08 

2013-06-19 23:35:55 

5.4 

200 

±14 

00032619013 

2013-06-27 20:17:08 

2013-06-27 23:54:56 

4.9 

170 

±13 

00032619014 

2013-07-03 00:03:17 

2013-07-03 22:18:55 

4.9 

165 

±13 

00032619015 

2013-07-10 09:26:47 

2013-07-10 14:40:55 

5.2 

327 

±18 

00032619016 

2013-07-17 11:14:55 

2013-07-17 21:16:54 

4.7 

299 

±17 

00032619017 

2013-07-24 11:46:07 

2013-07-24 20:00:55 

4.7 

166 

±13 

00032619018 

2013-07-31 06:55:12 

2013-07-31 18:35:55 

5.0 

241 

±16 

00032619019 

2013-08-07 00:48:43 

2013-08-07 22:04:55 

4.3 

192 

±14 

00032619020 

2013-08-14 15:33:40 

2013-08-14 19:12:55 

4.6 

207 

±15 


The source coordinates (J2000) are R.A. = 00**15™28?90 and 
Deck = —39°13H9”0 JStobbart et al.[[20()3i : the Galactic column 
densi ty along its direction is MjjGai = 1.37x10^^ cm“^ [Kalberla et alj 

I 2 OO. 5 I 1 . 


In this work, we present the results of a Swift monitoring 
campaign of the source NGC 55 ULXl, loca ted in a Magellanic 
barre d spiral galaxy at a distance of 1.78 Mpc ( [Karachentsev et aU 
[ 2 OO 3 I 1 and with X-ray luminosity in the range 8 x 10®® — 2 x 
10®® erg s“^. This low-luminosity ULX is peculiar because in 
two XMM-Newton observations it showed e nergy-dependent d ips, 
not commonly observed in ULXs (see also [Grise et al.[[2013l and 
[LinetalJlroIa for possible dips in NGC 5408 X-1 and 2XMM 
J125048.6-(-41074 3, respectively) . The dips lasted for a few hun¬ 
dreds of seconds ( [Stobbart et alJl2004) and appeared to be corre¬ 
lated with time variability at high energy. A later XMM-Newton ob¬ 
servation caught the source in a low flux state where dips were not 
present and the spectral properties were different (Pintore et al in 
prep.). The dips and the associated variability at high energy may be 
related to turbulences or blobs of optically thick matter in the wind 
that occasionally encounter our line of sight and mask the emitting 
inner regions (Middleton, Heil, Pintore, Walton and Roberts, sub¬ 
mitted). The study of this source appears then crucial to assess both 
the spectral variability and the onset of powerful outflows in ULXs. 
This fact, along with the goal of investigating the source long-term 
flux variability, motivated our request for a Swift monitoring cam¬ 
paign. Besides the Swift observations, we analyzed and present here 
also two previously unpublished Chandra observations. 


The paper is structured as follows: in Sectionl^we summarize 
the data selection and reduction procedures, in Sectionj^we present 
the results of our timing and spectral analysis of NGC 55 ULXl, 
and in Section|4]we discuss them. 


2 DATA REDUCTION 
2.1 Chandra observations 

We analysed the only two archival Chandra observations of 
NGC 55 ULXl (see Table[I), taken in 2001 September and 2004 
June wit h the Advanced CCD Imaging Spectrometer Imaging array 
fACIS-I: [Garmire et al.[[2003h . The CCDs were operated in Timed 
Exposure mode and in full frame (readout time: 3.2 s), and the 
Eaint telemetry mode was used. The data were reprocessed with 
the Chandra Interactive Analysis of Observat ions software pack¬ 
age (CIAO, version 4.6; [Eruscione et alJlrOOfih and the calibration 
files in the CALDB release 4.5.9. 

Despite the source NGC 55 ULXl was not at the aimpoint (4’ 
and 7’ off-axis in the first and second observation, respectively), in 
both observations the count rate was high enough to make pile-up 
effects in the ACIS detector not negligible. A ‘pile-up map’ cre¬ 
ated with the CIAO tool PILEUP_MAP confirmed the presence of 
pile-up. Owing to the sharp Chandra point-spread function (PSE), 
correcting for pile-up discarding the data in the core of the PSE 
would result in a stron g loss of cou nts. Therefore we used instead 
the PILE-UP model by Davis[ ( [200 ih (see sectionj^. The Chandra 
spectra should be treated carefully: since the pile-up strongly de¬ 
pends on count rate, the pile-up fraction was ~ 20% in the first ob¬ 
servation and ~ 40% in the second. We considered the pile-up level 
of the latter observation too high. The spectral mod ifications may 
be large enough that also the convolution model by [Davis[ ( [200111 
may not correctly describe them, making the spectral parameters 
not reliable. Therefore we splitted the second observation in two 
parts with the same duration and we spectrally re-analyzed the sec¬ 
ond part, when the source count rate was lower and pile-up was 
mitigated. The count rate was 0.34 cts s“^ and the pile-up grade 
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Figure 1. Light curve of the two Chandra observations (Obs.ID. 2255 and 4744 at the left and right, respectively) in the 0.3-10 keV energy band, sampled 
with a bin size of 100s. 



Figure 2. Swift light curve from the whole data set, with bin size of 4 days. 
We show the observations stacked in the Swift “group 1” and “group 2” 
spectra (see text) with blue and green colours, respectively. 


migration (a) decreased to 0.1 — 0.2, although these values are 
poorly constrained (see Table|2j- 

In observation 2255, the source counts were extracted from an 
elliptical region with semi-axes of 5” and 3”, while the background 
was estimated from an annulus with radii 8” and 15”. In observa¬ 
tion 4744, the semi-axes of the elliptical source region were 6” and 
9”, and the radii of the background annulus 15” and 25”. Back¬ 
ground regions were selected close to the source and we verified 
that the selected CCD regions were free of emitting X-ray sources. 


Swift. We collected approximately one 5-ks long observation per 
week in ‘Photon Counting’ mode (PC, t he XRT full imag ing mode, 
which has a readout time of 2.507 s; see lHill et alJl2004 ). 

The XRT data were processed with XRTPIPELINE (in the HEA- 
SOET software package version 6.14), filtered and screened with 
standard criteria. The source counts were extracted within a 20- 
pixel radius (~ 47”; one XRT pixel corresponds to ~2.36 arcsecQ), 
while the background was estimated from an annulus with internal 
(external) radius of 60 (100) pixels, i.e. 14l”(236”). We veri¬ 
fied that the background regions were in a CCD region free of other 
X-ray sources. For the spectroscopy, we used the standard CALDB 
spectral redistribution matrices, while the ancillary response files 
were generated with XRTMKARE, which accounts for different ex¬ 
traction regions, vignetting and point spread function corrections. 

The lightcurves were corrected using the XRTLCCORR tool, 
which takes into account Point Spread Function (PSF) corrections, 
vignetting and photon losses due to bad pixels/columns inside the 
extraction region. The latter corrections are based on a raw instru¬ 
ment map, constituted by a number of extensions corresponding 
to different time intervals, and its spatial dimension is arranged to 
cover the extraction region. 

The statistics of each individual Swift observation is poor and 
did not allow us to well constrain the spectral properties of NGC 
55 ULXl adopting the models described in the next sections. How¬ 
ever, the hardness ratio (1.5-10 keV/0.3-1.5 keV ) did not provide 
indications of statistically significant spectral variability between 
observations. Therefore, in order to improve the Signal-to-Noise 
ratio, we stacked all Swift observations according to their flux. We 
obtain two spectra (labeled as Swift “group 1” and “group 2”) com¬ 
prising 18 low flux and 2 high flux (Ohs. Id. 9009 and 9010) obser¬ 
vations, respectively (see the lightcurve in Figure 1). 


2.2 Swift observations 

We monitored NGC 55 ULXl be tween 2013 April an d August 

using the X-Ray Telescope (XRT: [Burrows et alJuOO^ on-board ^ http://swift.gsfc.nasa.gov/analysis/xrt_swguide_vl_2.pdf 
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Figure 3. A zoom of the light curve, in the 0.3-10 keV energy band, sampled with a bin size of 100 seconds in the first {left, center) and 50 seconds in the 
second (right) Chandra observation. There may be marginal hints of dips at the times marked with black arrows. 
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Figure 4. An example of three Swift observations in which marginal hints of dip episodes are seen (indicated by red arrows). The bin time of the lightcurves 
(in 0.3-10 keV energy band) is 300 seconds. 
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3 DATA ANALYSIS 
3.1 Timing analysis 

A peculiar property of NGC 55 UL Xl is the existence o f narrow 
dips in its XMM-Newton lightcurves. IStobbart et alj 12004) showed 
that, on average, the dips last about 100-300s (although a few dips 
appeared to be longer), with a depth at the level of ~ 80 — 90% of 
the persistent flux (~ 2.5 x 10“^^ erg cm~^ s“^) and more intense 
in the 2.0-4.5 keV energy band. We note that the average observed 
flux of NGC 55 ULXl in the Chandra and Swift observations is 
1.5 X 10“^^ erg cm“^ s“^, therefore lower than in the XMM- 
Newton data. If dips are characteristic of a specific high flux state of 
the source, we might not be able to detect any dip in the Chandra 
and Swift datasets. However, we can suppose that they are a ran¬ 
dom effect and present at different flux levels. In the XMM-Newton 
observations, we counted around 30 dips in a total exposure time 
of 51.9 ks. As they represent a small fraction of the total XMM- 
Newton exposure time, we calculated the expected number of dips 
detections Ndips in the Chandra and Swift observations according 
to the relation: Ndipa = (Wxmm/Ixmm) x t^^p, where Axmm is 
the total number of dips in the XMM-Newton observations, fxMM 


is the total XMM-Newton exposure time and t^xp is the Chandra 
or Swift exposure time. We therefore expected to detect 34 dips 
in the first Chandra observation (59.4 ks), ~ 5 dips in the second 
Chandra observation (9.6 ks), 2-3 dips in any single Swift snapshot 
(~ 4.5 ks) or 55 dips during the total Swift exposure time (95 
ks). 

The Chandra observations showed an increase in the count 
rate of about a factor of 2-3 between the two epochs (Fig. [TJ, and 
an average count rate in the 0.3-10 keV energy band of ~ 0.15 
cts“^ and ~ 0.4 cts s“^ for the first and second observation, re¬ 
spectively. Assuming to observe occultations similar to those of the 
XMM-Newton observations (depth at the level of 80 — 90% of 
the persistent flux), we should detect dips at ~ 0.03 and ~ 0.1 
cts s“^ in the first and second Chandra observation, respectively. 
We sampled the 0.3-10 keV lightcurve of the first Chandra ob¬ 
servation with 100s time bins as it assured the minimum binning to 
obtain Gaussian statistics in each time bin. We found only marginal 
indications of a limited number of dips, although the uncertainties 
due to the poor quality of the data are large (-^40% at the lowest 
rates; see some examples in Figure]^ left, center). The lightcurve 
of the second Chandra observation was instead sampled with 50s, 
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as it assured the minimum binning to obtain Gaussian statistics in 
each time bin. However, also in this case, we found only marginal 
evidence of dips (~60% error bars; see Figure]^ ngbf). 

The Swift observations generally showed variability of less 
than a factor of 2 (except for two high flux observations, that 
reached ~ 2.4 x erg cm“^ s“^; Fig.j^. Assuming a per¬ 

sistent mean rate between ~0.06 and ~0.1 cts s“^, we expected 
the depth of the dips to be at ~ 0.01 — 0.02 cts s“^, with a corre¬ 
sponding error bar of ~60-80%. The relatively low count rate de¬ 
tected by Swift did not allow us to investigate the existence of dips 
on timescales lower than 300s, as it represents the best compro¬ 
mise to reach gaussian statistics in almost each time bin and have 
small enough error bars. However, this implies that only the longest 
dips could be detected in these datasets. None of such deep dips at 
~ 0.02 — 0.01 cts s“^ was reliably detected, although we found 
possible hints in the lightcurve of a few observations where short 
drops in the count rate were observed (some examples are given in 
Fig.0. 

Dipping episodes in Galactic accreting sources are usually 
characterised by a spectral variability which can be clearly revealed 
by inspection of the hardness ratios. However, the large error bars 
of the hardness ratios in our Swift and Chandra observations do not 
allow us to clearly detect them and further investigate their prop¬ 
erties. Finally, we checked also the existenc e of dips in the 2.0- 
4.5 keV energy band (the same band in which l^obbart et ^|2004| 
found the strongest evidence) in both the Chandra and Swift obser¬ 
vations, but the statistics was too poor. 


3.2 Spectral analysis 


We analysed all the Chandra and Swift spec tra obtained a s de¬ 
scribed in Section musing XSPEC v. 12.6.0 ([a rnau 3 ll996h and 
grouping them to obtain at least 20 counts per bin. Two absorp¬ 
tion components (TBABS in XSPEC) were considered in all spec¬ 
tral fits: the first was introduced to model the Galactic absorp¬ 
tion al ong the direction of the source and hxed at 1.37 x 10^° 
cm~^ jKalberla et ahltZOtlSn . The second was allowed to vary and 
describes the absorption in the local environment of the source. 
The luminosities were estimated from the flux calculated from the 
CFLU X convolution model in X SPEC, assuming a distance of 1.78 
Mpc dKarachentsev et al.l2003h . 

_ E rom the analysis of the XMM-Newton data dGladstone et al.l 

120091: ISutton et al.f|2013l) , NGC 55 ULXl showed a complex 
continuum that can be described by a two-components model, 
consisting of a multicolour blackbody disc and a Comptonization 
component with a roll-over at high energy. However, this model 
was difficult to apply to the Chandra and the stacked Swift 
spectra because the data quality was poorer and did not allow 
us to constrain the spectral parameters. We then fitted the soft 
component with a single blackbody model (BBODY in XSPEC), 
neglecting that the emission of the outer disc/wind may come 
from layers with different ionisation and temperature. For similar 
reasons, the broad high energy spectrum, often described in 
terms of a Comptoniz a tion compon e nt (e. g 
[Gladstone et al.l 120091 : ISutton et alj l2013l : 


Stobbart et al 

20061 

Pintore et al.l 

2014t) 


was approximated with a multicolour blackbody disc (DISKBB in 
XSPEC). Thus, for both the Chandra and stacked S wift spectra we 
adopt ed, as common reference model, a DISKBB dMitsuda et al.l 
Il984h component plus an additional BBODY component (when 
needed). 


3.3 Swift spectra 

The “Group 2” spectrum does not need an additional BBODY com¬ 
ponent, as a single disc component, with a temperature of 0.6 keV 
provides a good fit to the data (see Table 1^. We note also that a 
DISKBB model is marginally favoured in comparison with a single 
POWERLAW component (x^/dof=44.39/34). Because of the disc¬ 
like spectral shape, we tried to fit the “Group 2” spectrum also with 
a modified disc model (DISKPBB. lMineshige et alj[i994l) . in which 
the radial dependence of the temperature goes as Tuisc oc r~^, and 
p is 0.75 for a standard disc and 0.5 for a slim disc. This model 
provided a good fit but its parameters turned out to be totally un¬ 
constrained. 

The best fit of the “Group 1” spectrum was instead obtained 
with a two-components model, a blackbody with a temperature of 
0.13 keV and a high energy DISKBB component with a tempera¬ 
ture of 0.56 keV (see Table^. As the mean flux estimated from 
the “Group 2” spectrum was about a factor of 4 lower than that 
of the “Group 1” Swift spectrum, we tested the possibility that the 
soft blackbody component might not be detected in the ’’Group 2” 
spectrum because of the poor counting statistics. We then added to 
the DISKBB component a BBODY model with fixed spectral param¬ 
eters, equal to those found in the “group 1” spectrum (see below) 
and letting the column density, the DISKBB temperature and nor¬ 
malization free to vary. We found that the additional component 
was acceptable and led only to a change in column density which 
grew from 0.13 atoms cm“^ to 0.49 atoms cm“^ (Table^. 

Comparing the spectral parameters of the two Swift spectra, 
we find that, while the total absorbed 0.3-10 keV luminosity 
increased (from 5 — 9 x 10^® erg s“^ to 1 x 10®® erg s“®) 
with the mean count rate of the stacked spectra, the unabsorbed 
luminosity did the opposite. The intrinsic source luminosity (i.e. 
unabsorbed) was about a factor of 4 higher in the “group 1” 
spectrum than in the “group 2” spectrum (~ 4 x 10®® erg s“® 
vs ~ 1 — 3 x 10®® erg s“®). This may be explained in terms 
of the higher column density (Ah) of the “group 1” spectrum 
(~ 0.6 x 10®® cm~®), about a factor of 6 higher than in the “group 
2” spectrum (~ 0.1 x 10®® atoms cm“®). However, comparing 
with the alternative two-components fit of the “group 2” spectrum, 
the luminosity of the “Group 1” spectrum is only a factor 1.5 
higher and the column density only 20% larger. 


3.4 Chandra spectra 

The parameters of the DISKBB and BBODY components of the 
first observation are kTub ~ 0.18 keV and kTdisc ~ 0.64 keV, 
slightly higher than those of the Swift spectra (but consistent with 
them within the errors). The spectral parameters of the second 
observation are less constrained. The column density is consistent 
with 0 within the errors, while the temperatures of the BBODY and 
DISKBB components (~ 0.15 and ~ 0.6 keV) are comparable 
to those of the first observation. We note that the intrinsic (i.e. 
unabsorbed) luminosity of the second Chandra observation is 
clearly the highest of the sample, an order of magnitude higher 
than the intrinsic luminosity of the “group 2” Swift spectrum. 

There is some hint of spectral variability between the Swift 
and Chandra observations, although most of it is accounted for 
by variations in column density (spanning the range between 
~1.5xl0®® atoms cm“® and 6x10®® atoms cm“®) and nor- 
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Table 2. Best fitting spectral parameters of NGC 55 ULX-1 obtained with an absorbed BBODY+DISKBB model. The error bars refer to the 90% confidence 
level. 



Nh“ 

kThb* 

kTdiBc'" 

/obs [0.3-10 keV]'' 

Lx [0.3-10 keV] ^ 

idi=kbb[0.3-10keV]/ 

X^/dof 


/'* 


10^^ atoms 

keV 

keV 

10" 

erg s~^ 

ergs-l 

10®*’ ergs-l 



% 

Chandra (Obs.ID. 02255) 

n 04+0-09 

Qg 

n 1 c+O-03 

n fi4+0-05 

U.D4_q Qg 


1 1 Q + 0.10 

. q + 0.10 
-'•■'^- 0.20 


141.44/149 

0 90+0-10 
U.2U_q 

0.85 

Chandra (Obs.ID. 04744)’ 


f) 1 c+O-OQ 

0 fil +0-08 
U-Di_o,io 


1 70 + 0.15 
^•'^- 0.20 

14.01^:2 


46.19/40 

0 1 O+0.23 
U.iO-Q 13 

0.85 

Swift, “Group 1" 

o.eolEi:;? 

0 1 Q + 0-02 

'J.lO_o.o2 

0.56iro4 


111 +0.03 

q Q + 1.30 

3.oi?:“ 

148.90/124 

- 

- 

Swift, “Group 2” 

f ) 1 O + 0.07 

u.iO-o.oe 

n 4Q+0-05 
U.4y_Q Q4 

0.13(/)- 

0.60l“:“ 

0.59l“:“ 


1 e /-)+ 0 . 2 O 

1-JU_q 20 

^•+U_Q 20 

1 1 + 0.20 
-'•■-'•- 0.10 
ry q + 0.20 
^■°- 0.20 

1 1 + 0-2 
-‘-■-‘-- 0.1 

1 q + 0.2 

-‘-■'^- 0.1 

39.45/34 

41.42/34 

: 

: 


^ Intrinsic column density of the neutral absorber; ^ Temperature of the blackbody component; ^ Inner disc temperature; ^ Observed flux in the 0.3-10 keV 
energy range; ® Unabsorbed luminosity in the 0.3-10 keV energy range, assuming a distance of 1.78 Mpc; ^ Unabsorbed disc luminosity in the 0.3-10 keV 
energy range, assuming a distance of 1.78 Mpc; Parameters of the JDPILEUP pile-up model as implemented in XSPEC: a is the grade-migration parameter 
(the probability that n events will be piled together but will still be retained after data filtering is / is the fraction of the PSF treated for pile-up and is 

required to be in the range 85-100%; 

* This observation has been split in two pails because of high pile-up; here we consider only the second pai1, where pile-up is lower. ** Alternative fit 
where the parameters of the BBODY component are fixed to those found in the “group I” Swift spectrum. 


Chandra, first obs. 


Chandra, second obs. 
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Figure 5. Top: unfolded spectra {E'^f{E)) of the first {top-left) and second {top-right) Chandra observations; bottom: unfolded spectra {E^ f{E)) of the 
stacked “group 1” {left) and “group 2” {right) Swift observations. The continuum is described by an absorbed blackbody (BBODY, dotted, red line) plus a 
multicolour disc (DISKBB, dashed, green line) or a simple DISKBB (see text). 
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kTbb (keV) 


Figure 6. Tentative linear correlation between the blackbody and disc tem¬ 
perature (taking into account the alternative fit of the “group 2” spectrum). 


malization of the spectral components. Using the alternative 
two-components fit of the “group 2” spectrum, we find marginal 
evidence of a linear correlation between the temperatures of 
the blackbody and disc components (kTjjisc = (1-4 ± 0.4) 
kT(,j, + (0.38 ± 0.06); see Fig. |^, with a regression coefficient 
of 0.92 and a Spearman coefficient of 0.8. This may suggest that 
the evolution of the two components is coupled, although it is not 
clear how it relates to changes in the mass accretion rate. 

Finally, in both the Swift and Chandra spectra, there is 
some indication (not statistically significant, though) of an ex¬ 
cess at ~ 1 keV in the residuals, that may be associated to 
an emission feature. Other ULXs showed a similar puzzling 
feature in their spectra, that has been interpreted as emission 
from a plasma close t o the source (i.e. StrohmaYgr_et_aL 2007|; 
Middleton_et_al2 20n|^ Caballero-Garcia et alj 20131 ; Sutton et^ 


20l3t 


Pintore et al 


I Middleton et alj|2(^ a). 


20 ^ or mismodeling of absorption features 


4 DISCUSSION AND CONCLUSIONS 

We analysed two Chandra ACIS-I observations and 20 Swift obser¬ 
vations, taken in PC mode, of the Ultraluminous X-ray source NGC 
55 ULXl. We investigated the spectral variability of the source, 
stacking the Swift observations in two spectra according to the ob¬ 
served count rate, with the aim to improve the signal-to-noise ra¬ 
tio. For both the stacked Swift spectra and the Chandra spectra we 
adopted a model consisting of a blackbody (for the soft compo¬ 
nent) and a multicolour accretion disc (for the hard component). 


lar to that observed in other ULXs (e.g. 

Kaiava & Poutanen 1 20091; 

iFeng & Kaaretll2009l; ISutton et al.1 12013 

; Pintore et al. 20l4). We 


can account for this spectral variability mostly in terms of changes 
in both normalization and column density of the two model com¬ 
ponents. 

Despite detailed timing information are missing, as many 
other ULXs, NGC 55 ULXl escapes an interpretation in 
terms of transitions between canonical sub-Eddington states as 
those observed in accreting Galactic BIT X-ray binaries (e.g. 
iMcClintock & Remillardll2006l). Althoug h different scenarios have 


been proposed (e.g. Miller et al.l l20l4) . many of the observed 


properties can be qualitatively accounted for by the expected 
behaviour of an accretion disc accreting above Eddington (e.g. 
Gladstone et al.ll2009l; IZam pieri & Roberts! l20()9l ; iMiddleton et al.l 


2 OI 2 I; ISutton et al.l |2013|; ' Ipintore et al.l ' |2014 . In such condi¬ 


tions, advection of energy takes place in the accretion disc, 
the structure of which is modified with respect to that of 
a standard disc jAbram owicz et akl Il988l ; IWatarai et alj 1200 ll ; 
lOhsuga & Mineshig3 2011 1. In addition, part of the accretion en¬ 
ergy is removed in forms of ra diative outflows that originate in the 
inner regions of the disc (e.g. [Poutanen et al.ll200^ ; lOhsuga et al.1 
l2009h . Their opening angle (measured from the rotation axis 
of the disc) decrea ses when the accretion rate increases (e.g. 
iTakeuchi et al.ll2oI4 . The density of these outflows is sufficiently 
high that they are optically thick and produce thermal emission. 
The characteristic temperature of the photosphere in such winds is 
consistent with that of the soft component often detected in ULX 
spectra. If the angle between our line of sight and the rotation 
axis of the accretion disc is sufficiently large, short term variabil¬ 
ity jHeil et al.ll2009ll can be induced by o ptically thick clouds of 
matter that interse ct the edge of the winds ^Middleton et alJlrOl ll ; 
ISutton et al.ll2013l) . In case the wind is sufficiently dense and ex¬ 
tended, its outer neutral phase can also partially obscure the high 
energy emission from the inner disc and its own soft emission, 
leading t o an increase in th e local absorption and column den¬ 
sity (e.g. IPintore et al.ll2014h . Narrow absorption lines from this 
warm phase of the wind are still not yet convincingly detected 
dWalton et al.ll20l3) . although broad, blue-shifted X-ray emission 
features of ionised spe cies of abundant elements may be expected 
dMiddleton et alj2014h a. In the following, we then interpret the ob¬ 
served spectral evolution of NGC 55 ULXl as induced by varia¬ 
tions in accretion rate, outflow rate and intrinsic absorption assum¬ 
ing that the source is in a super-Eddington accretion state of this 
type and that the inclination angle to the source is large. 


We then start from the second Chandra observation, the in¬ 
trinsic luminosity of which reaches 10“^° erg s“^ (Table|^, a value 
higher than that of any XMM-Newton observation. Both the soft 
and hard emission are stronger than in any previous observation. 
We then suggest that, at the time of this observation, the accretion 
rate is the highest, and that the outflow is extremely powerful with 
a very extended photosphere. In the first Chandra observation the 
spectrum shows still evidence of a powerful soft component, but 
with a total intrinsic X-ray luminosity of ~ 1.3 x 10®® erg s“^. 
The accretion rate is then probably lower than in the second Chan¬ 
dra observation. The “group 1” Swift observations represent a sort 
of intermediate state in which the accretion rate is probably smaller 
than that of the second Chandra observation, but higher than that in 
the first. While the outflow was still strong, the ionising flux from 
the disc decreased and then the fraction of neutral wind material 
became higher, absorbing more effectively its own emission and 
that of the inner regions. We note that the column density of these 
observations is ~ 0.6 x 10®^ erg s“^, the highest amongst all the 
spectra considered here. The absorbed luminosity (and count rate) 
is so low because the view of the central regions is highly sup¬ 
pressed by this obscuring low-ionisation component. Finally, the 
interpretation of the state of the source during the “group 2” Swift 
observations depends somewhat on the spectral fit considered (see 
Table|3- Assuming the single component disc fit, the intrinsic low 
luminosity and intrinsic column density, on one side, and the lack 
of the soft component on the other, may suggest that the accretion 
rate is at the lowest level. According to nu merical simulations (e.g. 
[Poutanen et alj|200^ ; lOhsuga et alJl200^ . this would represent a 
state in which the wind photosphere is not extended and the wind 
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opening angle decreases to the point that it is no longer intersecting 
our line of sight towards the central regions. Considering the two- 
components fit, the situation is less clear, hut may be suggestive of 
a relatively weak wind with a still well developed neutral fraction. 
As the two “group 2” Swift observations were collected within 20 
days after the first low flux interval (see FigO, this second sce¬ 
nario may be preferable unless we allow for a total suppression of 
the wind in such a small timescale. 

One of the tests of the super-Eddington scenario is the de¬ 
tection of absorption features originating in the warm wind. The 
counting statistics of both the Swift and Chandra data did not al¬ 
low us to perform a detailed analysis, but we note that there are 
marginally sig nificant residuals in the spectral fits at ~ 1 keV in 
all the spectra. [Middleton et al.l J^Oldh a interpreted features of this 
type as arising in the wind. Deeper observations and higher quality 
spectra will help to confirm these claims. 

Finally, we comment on the dips in NGC 55 ULXl. One 
of the aims of the Swift monitoring was to unveil possible flux- 
dependent dips. In the XMM-Newton data, dips lasting for 100- 
300s were observed only during the high flux states, and they 
were explained as caused by optically thick clouds of matter in the 
sourc e environment, that f rom time to time encounter our line of 
sight dStobbart et alj|2003l . If the mechanism for producing them 
is similar to that of the Low Mass X-ray binary systems (LMXBs), 
then the source i s seen at relatively high inclination (> 65°, e.g. 
[White et al.li^ . However, although we have observed the source 
spanning a rather large interval in flux (a factor of 3-4 in both the 
Chandra and Swift data), no clear dip episodes were observed but 
only marginal drops in the flux. Assuming occultation of 80-90% 
of the persistent flux, for the given exposure times of the Chandra 
and Swift observations, we would expect to detect a large number 
of dips. However the counting statistics of these datasets was poor, 
with large uncertainties (~40-60%) at the flux level of the candi¬ 
date dips: for the Swift data, we could investigate only dips longer 
than 300s, while in the Chandra data we could investigate shorter 
timescales. In some Swift observations and in the two Chandra ob¬ 
servations, we found marginal hints of narrow drops in the count 
rate although most of them are consistent with the mean within 
3 (t. A n analysis of the lightc urves in the 2.0-4.5 keV energy band, 
where [^obbart et ^ i2004h found the strongest evidence of dips, 
showed that the counting statistics was, even more, too poor to al¬ 
low detection of dips. In addition, the hardness ratios do not show 
any clear spectral change between dips and ‘persistent’ episodes. 

The lack of evidence of detectable dips in our data can there¬ 
fore be attributed either to poor counting statistics or to the ran¬ 
dom occurrence of favourable/unfavourable obscuring conditions 
by blobs in the wind. The latter is ultimately related to the un¬ 
known physical interp lay between the accr etion phase and the 
disc/wind propagation ([Middleton et al.|[20l^ ) in NGC 55 ULXl. 
Indeed, also the spectral variability of the soft component suggests 
that the wind, and its ionisation, could be highly variable on time 
scales from a few ks (the duration of the dips in the XMM-Newton 
lightcurve) to several days (high flux Swift observations). 

Further spectral investigations, complemented with a timing 
analysis, using high quality observations will allow us to further 
improve our understanding of NGC 55 ULXl and the role of out¬ 
flows in ULXs. 
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